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We explore the phase diagram of interacting spin-1/2 systems in the presence of anisotropic inter¬ 
actions, spontaneous decay and driving. We hnd a rich phase diagram featuring a limit cycle phase 
in which the magnetization oscillates in time. We analyze the spatio-temporal fluctuations of this 
limit cycle phase at the Gaussian level, and show that spatial fluctuations lead to quasi-long-range 
limit cycle ordering for dimension d — 2. This result can be interpreted in terms of a spatio-temporal 
Goldstone mode corresponding to phase fluctuations of the limit cycle. We also demonstrate that the 
limit-cycle phase exhibits an asymmetric power spectrum measurable in fluorescence experiments. 
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A quantum system that is coherently driven and con¬ 
nected to a heat bath eventually reaches a steady state; 
when the system is macroscopic, this steady state can be 
ordered in the sense of spontaneously breaking a symme¬ 
try mm- The patterns of steady-state ordering are, in 
general, qualitatively different from those of the equilib¬ 
rium phase diagrams: for instance, the interplay between 
coherent and dissipative dynamics can stabilize staggered 
phases that are absent in equilibrium [3], and by engi¬ 
neering the dissipative terms one can optically pump a 
many-body system into a pure state mm- The forms of 
steady-state order [MIS] that have been investigated to 
date mostly fit into equilibrium paradigms, such as the 
Landau paradigm of spontaneous symmetry breaking or 
the paradigm of topological order. However, the depar¬ 
ture from equilibrium allows one to realize novel types of 
order—specifically, limit cycles (LG) pTJ[2Q] which break 
time-translation invariance—that have no obvious equi¬ 
librium counterpart. Although previous works have pre¬ 
dicted a LG phase within mean-field theory p^H2Q] . it is 
an open question whether such a phase really exists or is 
an artifact of mean-field theory. 

In this work, we show that the LG phase exists with 
long-range order in three and higher dimensions and 
quasi-long-range order in two dimensions. We study 
a paradigm model of interacting spins—the anisotropic 
spin-1/2 Heisenberg model in a transverse field—and find 
a regime where it exhibits a limit-cycle phase. We dis¬ 
cuss the origin of the LG and study the effects of quantum 
fluctuations to go beyond previous mean-field works. We 
show by explicit, microscopic calculation that the spon¬ 
taneous breaking of the continuous time-translation sym¬ 
metry is reflected in the presence of a gapless “Goldstone” 
mode, and that this gapless mode prevents global time- 
translation-symmetry-breaking in one or two dimensions. 
We then discuss this LG at a more phenomenological 
level, noting its unusual implications, e.g., that the tem¬ 
poral ordering of the LG phase gives rise to an asymmet¬ 
ric dynamical power spectrum of emitted photons. These 
predictions are straightforward to test in experiments 


with trapped ions [2TH23] or Rydberg atoms [3 HI HU- 

We want to contrast our paper with recent works 
on time crystals which also spontaneously break time- 
translational symmetry [261431] . A LG is the nonequilib¬ 
rium steady state of a master equation, whereas a time 
crystal is the ground state of a Hamiltonian; thus, LGs 
are unaffected by the no-go theorems concerning equilib¬ 
rium time crystals [30] . 

Model and phase diagram .—Gonsider a d-dimensional 
square lattice of spins with nearest-neighbor exchange in¬ 
teraction. Each spin is resonantly driven by a laser and 
experiences spontaneous decay that incoherently flips the 
spin from up to down. The resultant many-body dynam¬ 
ics is governed by a master equation for the system’s 



FIG. 1: (Color online) (a) Mean-held phase diagram of the 
driven-dissipative spin system for {Jy,Jz) /y = (6,2), high¬ 
lighting the limit cycle and its vicinity. The nonequilib¬ 
rium spin system has different phases: uniform (Ui), bistable 
(U 1 /U 2 ), staggered (S), and LC phases. LC/U 2 denotes phase 
coexistence. In the LC phase, the spin system undergoes 
self-sustaining collective oscillations in time. The solid and 
dashed lines represent continuous and discontinuous phase 
transitions, respectively, (b) An example trajectory of the 
time-dependent LC at (Ut, Jx, Jy, Jz) h — (I?“f?65 2) that 
breaks both the time-translational and sublattice symmetry. 
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FIG. 2 : Zoomed-out phase diagram showing different phases 
for (Jj^, J^)/ 7 = ( 6 , 2 ). 


density matrix p: 


p = -i[Hj+HQ,p]-^'^ 

i 


2 Wt ^P]- P^t 


( 1 ) 




where erf is the a component Pauli matrix for spin i. Hj 
and describe the Heisenberg interaction and coher¬ 
ent driving, respectively. 7 is the spontaneous decay rate, 
experimentally introduced via optical pumping. The ex¬ 
change interaction can be realized in trapped ions via 
motional sidebands [2TH23] . and Rydberg systems based 
on optical adiabatic elimination HlEilEl]. J and 7 can 
be engineered to be on the order of 1 kHz with trapped 
ions or 100 kHz with Rydberg atoms. In the following, 
we focus on the parameter regime J ^ ^ 7 . 

The resulting dynamics can be understood in terms 
of the interplay between spontaneous decay, drive, and 
spin interaction. When there is no interaction, a spin is 
driven around the x axis and equilibrates in the lower 
yz plane (so that (a^) = 0, (a^) 7 ^ 0, (a^) < 0). In the 
presence of interaction, each spin precesses about an ef¬ 
fective field neighboring spins. 

This effective field is established self-consistently, leading 
to the possibility of ordered states. When the coupling 
is anisotropic ^ ~Jy)i fhe effective field becomes 
almost perpendicular to the spin direction, and the pre¬ 
cession effect is stronger. Therefore, we expect a richer 
phase diagram in the anisotropic coupling regime. Note 
that this system does not possess any Z 2 symmetry (like 
erf^a^ —erf,—erf) due to the presence of the drive. 

The only symmetries of the master equation are a con¬ 
tinuous symmetry under time translation and a discrete 
symmetry under spatial (lattice) translation. We restrict 
our analysis to bipartite lattices; thus, the mean field on 
sublattice A is due to the magnetization on sublattice B 
and vice versa. 

A typical scan of the steady-state phase diagram is 
shown in Figs, [^a) and|^ The calculation is based on 


FIG. 3: (Golor online) (a) Phase transitions of the LG phase 
as a function Jx at {Q, Jy, Jz)/^ = (0.6, 6 , 2). The sublattice 
spin order parameter develops bistability, Hopf bifurcation 
and sublattice order. The LG phase has a continuous and 
discontinuous phase transition to the Ui and S phase, respec¬ 
tively. Vertical bars represent the amplitudes of the sublattice 
spin oscillations. Gircles and solid squares denote the two sub¬ 
lattice stationary states. The LG amplitude shows a scaling of 
near the continuous LG-to-Ui transition, (b) LG 
period T as we vary H for fixed (Jx, Jy, Jz)H — (~6.37, 6 , 2). 
The divergence of period scales as T ~ (( 5 H/ 7 )~^^^ as we 
approach the tricritical point at (Jcc,c, f^c )/7 ~ (—6.37,0.37). 


sublattice mean-field theory and linear stability analy¬ 
sis to obtain the stable fixed points of Eq. Q . There 
are (i) three qualitatively different uniform steady states 
(Ui, 2 , 3 ) that do not break any symmetry and predomi¬ 
nate at large drive and/or isotropic couplings; (ii) a stag¬ 
gered phase (S), which is time-independent but breaks 
the sublattice symmetry, and dominates for strongly 
anisotropic coupling; and (iii) a LC phase, which breaks 
both sublattice and time-translation symmetries. 

The uniform phases are related to the paramagnetic 
(U 3 ) and the two ferromagnetic (Ui^ 2 ) steady states of 
the undriven system; in the presence of the drive, Z 2 sym¬ 
metry is explicitly broken so these phases are no longer 
distinct in terms of symmetry. The staggered phase is 
a descendant of the antiferromagnet; since it breaks lat¬ 
tice translation symmetry, it is quite different from the 
uniform phases; thus the uniform-staggered transition is 
second-order. Finally, the LC phase, which has no ana¬ 
log in the undriven system, emerges near the uniform- 
staggered phase boundary in the presence of a drive. It 
is intuitive that the LC arises in this regime, since even 
in the undriven system the relative orientation between 
the two sublattices is “soft” (i.e., highly susceptible) near 
the transition between the ferromagnet (Ui) and the an¬ 
tiferromagnet (S). 

Limit cycle: mean-field analysis .—In the LC phase, 
the magnetizations on the two sublattices oscillate with 
a relative phase of tt [Fig. [^b)]. Thus, the LC breaks 
both spatial and time translation symmetry. 

Figure]^ (a) illustrates the properties of the LC phase 
transitions. The plot shows the mean-field amplitude of 
Sx _ vary with fixed ^2, Jy and Jz- As 
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Jx decreases, the spin system goes from a bistable phase 
(U 1 /U 2 ) to a uniform (Ui) phase, and then exhibits a 
supercritical Hopf bifurcation [HEg to the LC phase at 
the critical point of Jx^c ~ — 6 .O 87 . A further decrease 
of Jx to ~ — 6.787 renders a first-order transition from 
LC phase to the staggered phase. Note that in this LC 
phase, there are no stable fixed points. The amplitude of 
the LC scales as near the second-order phase 

transition from LC to Ui, with 5Jx = Jx,c — Jx [33] . 

The intersection of continuous LC-Ui and discontinu¬ 
ous LC-S phase boundaries defines a LC tricritical point 
[Fig .[^a)]. Interestingly, the LC period diverges near this 
tricritical point [Fig. [^b)]. The reason stems from the 
imaginary part of stability eigenvalues that determine the 
frequency of LC [33]. Along the LC-Ui boundary, the two 
stability eigenvalues are purely imaginary complex con¬ 
jugates, corresponding to the Hopf bifurcation, while on 
the S-Ui boundary, they both vanish due to the sublat¬ 
tice symmetry. Therefore, the LC near their intersection 
has a divergent period [34] . 

Gaussian fluctuations. —We now go beyond mean- 
field theory and explore Gaussian fluctuations about the 
mean-field states. Because the mean-field LC spon¬ 
taneously breaks time-translation invariance, fluctua¬ 
tions about it obey an inhomogeneous Floquet equation 
[35l [36] . The correlation functions can be expressed as 
(see Supplemental Material 137]): 

C(/c, t) = A(/c, t)C{k, t) + C(/c, t)A^{k^ t) + D{k^ t), ( 2 ) 

where C{kfl) is the Fourier transform of the spatial 
correlation matrix t) = (crg( 0 , t)cr^/(r, t)) — 

(crg(0, t))(cr^/(r, t)) with L and L' denoting sublattices. 
A and D are the drift and the diffusion matrices, re¬ 
spectively [38|, and are both periodic in time due to de¬ 
pendence on the mean-field variable SL{t). So, C is also 
time-dependent. In contrast, for a stationary phase, the 
steady state solution for C does not have any time de¬ 
pendence. 

In the case of a LC, the correlation function has inter¬ 
esting dynamical behavior due to the time periodicity of 
A and D. Figure |^a) presents a typical plot of the mo¬ 
mentum and time dependence of the solution of Eq. <§ in 
the LC phase. Figure |^b) shows that the oscillating cor¬ 
relation function reaches a steady value for finite k, but 
diverges linearly in time at /c = 0. In the long time limit, 
the stroboscopic correlation goes as l/k‘^ for small mo¬ 
mentum [Fig. [^c)]. This infrared divergence precludes 
LC ordering for dimensions d < 2, and—as we discuss 
below—can be traced to the Goldstone mode [T] gener¬ 
ated by spontaneous time-translational symmetry break¬ 
ing in the LC phase. In 2 < d < 4, a similar analysis can 
be used to compute correlation lengths [Fig. [^d)], and 
to estimate the size of the fluctuation-dominated region 
around the onset of a LC via a Ginzburg criterion (see 
Supplemental Material m)- In three dimensions, the 


(a) C'l^(k,t) (b) 



Cl\{k = ko,t) 



0.2 


(c) 



FIG. 4: (Color online) (a-c) Dynamical spin correlations when 
(H, dec, Jy, Jz)ll — (0-^? —6.09,6, 2 ). (a) A typical plot of the 
sublattice correlation function as a function of momentum 
and time in the long-wavelength limit, (b) When /c = 0, 
the correlation diverges linearly with time, while it saturates 
to finite values for /c 7 ^ 0. The oscillations have the same 
period as the LC. (c) The stroboscopic correlation in the long 
time limit shows a 1/k^ dependence, reflecting the gapless 
excitation of the broken time-translational symmetry of the 
LC phase, (d) Correlation length ^ scales as ^ ~ 
near the continuous LC-Ui phase transition. 


critical region has a parameter-space width 5Jx 10 “^ 7 , 
which is of the order of 10 kHz in atomic experiments. 

Phenomenological description .—The analysis above 
has the advantage of being fully microscopic, thus en¬ 
abling comparisons with experiment; however, it does not 
give a transparent picture of the origin of the Goldstone 
mode. In the spirit of Landau theory, we now consider a 
general phenomenological description of a LG, known as 
its normal form [33| . The idea is to expand the order pa¬ 
rameter around the time-independent component Sc near 
the onset of a Hopf bifurcation: S{t) — Sc = A(t)p + c.c. 
with the vector p lying on the plane of oscillation. Fol¬ 
lowing standard procedures [33] , the mean-field equation 
can be rewritten near the LG-Ui boundary as: 


(iA(x, t) 
dt 


^A+ ^ + A{x,t) 

A{x,t) + C(t), (3) 


where A, T are the amplitude and period of the LG, re¬ 
spectively. We have augmented the normal form with 
a stiffness S against spatial fluctuations, as well as a 
Langevin noise term [39]. The stability of the LG re¬ 
quires the coefficients A G 5R (being linear in and /i 
to satisfy A Re /i < 0. Importantly, Eq. reveals the 
continuous U{1) symmetry {A Ae^'^) of the LG. 
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The normal form Eq. (§ can be used to derive the 
lower critical dimension by the following procedure (see 
Supplemental Material m- working in the rotating 
frame (^4 ^ expanding A in terms of phase 

fluctuations to obtain a phase-only action, and using 
this action to compute the correlation function of the 
phase of the LC. In momentum space, this correlator is 
(^k^-k) ^ 1/^^: and thus its inverse Fourier transform 
in 2D is 

(6>(r)6>(0)) ^ 7 'logr ^ {A{r)A{0))r^l/r^, (4) 

i.e. algebraic long-range order with 7 ' ^ ^/{SAq) and Aq 
being the LC amplitude. This corresponds to effectively 
thermal behavior with an effective temperature set by 
the noise term ( (which is proportional to 7 ). 

Photon temporal correlation .—The LC phase possesses 
novel dynamical features in the temporal fluctuations of 
spins, which can be seen in the emitted fluorescence. 
Each spin is coupled to the electromagnetic vacuum such 
that the time correlation of the fluorescence signal re¬ 
flects that of the spin. For a conventional stationary 
state, like the Ui phase, the two-time correlation function 
(cr+(0)cr“(r)) has a damped oscillation, corresponding to 
a power spectrum with a three-peak structure [Fig. da)]. 
This is because the spins are collectively driven by an ef¬ 
fective field and the dynamics are similar to those of a 
driven two-level system, so the photon power spectrum 
resembles a Mollow triplet [40] . 

However, the temporal dynamics are very different in 
the LC phase. A power spectrum is usually defined for 
a stationary state. Extending it to a time-dependent 
phase, the dynamical power spectrum is: 

T 

SL{t,uj) = —Re dr((7(^)o-i(i + r)))e*“'^, (5) 

Jo 



h 

0.3 

j 


FIG. 5: (Color online) Power spectra as a function of detun¬ 
ing A with respect to the spin natural frequency in stationary 
and LC phases, (a) At {Q, Jx, Jy, Jz)/"f = (0.6,—6, 6, 2), the 
system is in the time-independent Ui phase and the corre¬ 
sponding power spectrum has a standard Mollow triplet struc¬ 
ture. (b) The sublattice power spectrum for the LC phase at 
(D, Jx, Jy, Jz)h — (I? ^ ^5 2) and a particular time t. The 

asymmetry in A is caused by the complex-valued two-time 
correlation function in the LC phase. Both sublattice spectra 
have a time-dependence and they do not cancel each other, 
(c) The total dynamical power spectrum SA{t,(ju) + SB{t,(ju) 
oscillates in time with the same periodicity of the LC. 



To probe the dynamical power spectrum, one can per¬ 
form a time-resolved photon detection. The photon de¬ 
tection for a time T is proportional to 001 : 

1 

P{T,u)) oc -Re dt , ( 6 ) 


where the fluctuation is ((O 1 O 2 )) = {O 1 O 2 ) — ( 0 i)( 02 ), 
Jo is the conventional power spectrum coefficient m and 
L denotes sublattice. For a stationary phase, only the 
time difference r between the spin operators matters and 
the spectrum is independent of t. However, in the LC 
phase, {t)a^ {tr)) 7 ^ (cr^( 0 )cr^~(r)) and the spec¬ 
trum becomes periodic in t. To calculate it, we focus on 
the parameter regime far away from the phase boundary 
so that spatial fluctuations are negligible and mean-field 
theory is valid (see Supplemental Material 137]). 

Figurej^b) shows the power spectrum for the LC phase 
at a given time t after the system reaches a steady state. 
The spectrum is asymmetric in the detuning because 
(cr+(t)cr“(t + r)) is no longer real due to the loss of 
time-translational invariance. Moreover, as illustrated in 
Fig. I^c), the power spectrum is periodic in time due to 
the oscillations. These dynamical properties are in sharp 
contrast to fluorescence signals in a stationary system 
[Fig- da)]. 


where the average value is subtracted. A T-resolved 
photo-detection permits the measurement of (T + 
5T)P{T P 5T,uj) - TP{T,uj) ^6TxS{T,uj), when 6T is 
much shorter than the LC period. 

In the case of quasi-long-range order or short-range 
order with a large correlation length, the power spec¬ 
trum will become symmetric when averaged over long 
times. However, one expects that higher-order correla¬ 
tion functions, such as {{crL{t)crL{t + T)aL{t')aL{t' + r))) 
(with T ^ T and \t — t'\ ^ T) will retain signatures of 
limit-cycle behavior, i.e., they will exhibit oscillations at 
a period T. 

Conclusion. —Driven-dissipative spin systems with 
anisotropic interactions have a rich nonequilibrium phase 
diagram; the most remarkable feature of this phase di¬ 
agram is the presence of a limit-cycle phase. Although 
the limit cycle is inherently nonequilibrium (and thus ex¬ 
hibits an asymmetric power spectrum), this phase and its 
transitions can be understood in terms of familiar equi- 
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librium concepts such as correlation lengths, Goldstone 
modes, and critical dimensions. We believe this work 
opens up a fruitful new interface between the condensed 
matter and nonlinear dynamics communities. Possible 
future directions include studying the criticality of spa- 
tiotemporal oscillation patterns im, non-Markovian en¬ 
vironments [42], and non-Hermitian spin systems [43] . 

This work was supported by the NSF through a grant 
to ITAMP. C.K.C. acknowledges the support from the 
Croucher Foundation. 
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Supplemental Material: 

Limit cycle phase in driven-dissipative spin systems 

Ching-Kit Chan, Tony E. Lee, and Sarang Gopalakrishnan 


This supplement details the calculations of Gaussian theory and Floquet analysis for spatial correlations, a phe¬ 
nomenological derivation of the lower critical dimension and the photon temporal correlations in the LG phase. 


EQUATIONS OF MOTION AND MEAN-FIELD THEORY 

According to the master equation [Eq. Q] in the main text, the exact equations of motion for averages of one and 
two spin operators can be systematically written as: 

/X 

- (r“ + T^) (a>^) - 7 

/X 

where the subscript (n^m) and the superscript ((a,;d, 7 , (5) denote the site and the spin index, respectively. Here, 
= 7 and /i = The second equation only holds for n ^ m. We have also used the 

Levi-Givita symbol to simplify the notations. Both 5^^ and 8 nm are Kronecker deltas. 

The sublattice mean-field decoupling of Eq. Q (i.e. (crf+zx^n) ~ ~ ^b^a) results in a 6 -components 

self-consistent equation. A linear stability analysis of the steady state solution leads to the sublattice mean-field phase 
diagrams Fig.j^a) and Fig. in the main text. Note that there can be multisite instability such as spin-density-wave 
[ 1 ] that goes beyond the sublattice mean-field theory. We have checked that and the spin-density-wave instability 
occurs within the U 3 phase and does not affect the LG. 


(7) 

( 8 ) 


GAUSSIAN SPATIAL FLUCTUATION 

Within the Gaussian framework, we define the two-point correlation functions and approximate higher order cor¬ 
relation functions as: 


^3 \ -a -t 

T rr^ ) = rr rr^ 


^a,l3 

m 




n^m^l 




(7j T ^n,m ’ 


when n ^ m ^ 1 . 


(9) 

( 10 ) 


After a further introduce of sublattices L = ±1 (or A and B), the mean-field variable and the correlation function 
become = dg and Applying these approximations, Eq. 0 becomes: 


C 


L,L',r 


+ J E {^iCliL>,r+, + ^lLCZ^,r + <5-L,L'<5,+^,o7 - 7,4 - C27,o) } 

+ J E + ^-L'C2i',r + S-L,L>Sr+,,o^i, - Clio) }(,11) 


which is valid for r 7 ^ 0. This equation can be understood in the followings. The first line comes from the driving 
and the decay contributions. In the second and third lines, the first two terms come from the Gaussian decoupling, 
and the third term accounts for situations with equal site indices. For example, when r = —/x, terms like 
only consist of spins from two different sites and shall not be decoupled. In this case, should be replaced 

by + S^^]) instead. 
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be expressed as: 


To solve Eq. (11), we go to the momentum space. Defining ^ the Fourier transform of Eq. (11) can 


C q;,/3 _ I ^a,6 aP,5 , ^a,/3 

L L' k — L" k^L" L' k ' h" k^L' L" k L L' A;’ 


where A and D have the meanings of the drift and the diffusion matrix, respectively, with 


D 


a,^ _ 

L,L',k ~ 


5- 


L,L> 


d 


- CT2,af,) + {ie^^^dl + <5“^' - dial) 


( 12 ) 

(13) 

(14) 


We have dropped qualitatively unimportant terms ^ 0{l/d?) in the D expression. For a stationary phase, the steady 
state correlation can be obtained by putting ^ ~ However, in the LC phase, dg(t), A{t) and D(t) are 

time-periodic, and one has to solve the 36-components Eq. (12) numerically. 


FLOQUET ANALYSIS FOR THE SPATIAL FLUCTUATION 


We now extract the analytical properties of the solution of Eq. (12), which is an inhomogeneous linear differential 
equation with time-periodic coefficients. It can be rewritten as: 


Ck{t) = Mk{t) • Ck{t) -t- D/e(t), 


(15) 


where the vectors Ck{t) is the 36-components correlation functions, Mk{t) and Dk{t) contains the drift and the 
diffusion coefficients, respectively. Let T be the period of M(t)^ and the LC. According to the inhomogeneous 
Floquet theory m, the solution of Eq. (iT^ can be expressed as follows: 


Ck(t) = Xk{t) 


Ck{0) 


f 


dTX^\T)Dk{T) 


(16) 


where X^it) is known as the principal fundamental matrix of Eq. (15). The Floquet theorem says that the fundamental 
matrix has the following form: 


Xk{t)=Pk{t)e^’^\ 


(17) 


where Pk{t) is some T-period matrix and F is some time-independent matrix. Note that Xk{t) is obtained from the 
homogeneous equation Xk{t) = • Xj^{t), with an initial condition X/e(0) = I. 

The stability of Eq. (15) is governed by the eigenvalues of Xk{t = T). These eigenvalues, denoted by are 

known as the Floquet multipliers and jii are referred as the Floquet exponents. When jii < 0, the result is stable; 
while if jj^i > 0, the solution will grow exponentially in time. Figure [^presents the largest two Floquet exponents as 
a function of k in the LC phase. In the long-wavelength limit, ^ —k^ and /i 2 ^ “(^o + Physically, and 
fi 2 refer to the transverse and the longitudinal fluctuations, respectively, and ko ^ Therefore, for finite k values, 
the correlation function is stable and oscillatory in time [Fig. [^d)]. The k = 0 situation is special since the largest 
Floquet exponent is zero, corresponding to a Floquet instability. 


Correlation length, upper critical dimension and Ginzburg criterion 


With the help of Floquet exponents, we can obtain the analytical behavior of the correlation functions. One can 
show that for t = nT^ Eq. (16) can be rewritten as [3]: 


Ck{t = nT) = Xk{T)^Ck{0) + [Xk{T)+Xk{T)^ + ... +Xfc(T)"] TdrX^i(r)A(r). 

^0 


(18) 


For long time (n ^ 1), the above expression is dominated by the second inhomogeneous contribution. Diagonalizing 
Xk{T) = Uk{T)Kk{T)U^^{T) with {Ak{T))- = we can write in terms of the Floquet exponents: 


Cu{t = nT) =[Uu{T)], 




(n-|-l)T \ 
1 - eA'.T’ ) 


\T) f 
Jo 


dTXl\T)Dk{T) 


3 


(19) 




















(a) 

-/ii(k) 
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(b) 

-Ai2(k) 



FIG. 6: (Color online) The first and second largest Floquet exponents fii and for the correlation functions C{k,t) in the LC 
phase at (Q, Jy, Jz)l'y = (0.6, 6, 2). fii oc —k^ and 112 ~ —(^^ + ^ 0 ) in fhe small k limit, is the inverse correlation length 
which diverges as Jx/'J approaches the continuous phase transition. 


Using the k dependence of the dominant Floquet exponents, each stroboscopic correlation function in the long- 
wavelength limit can be expressed in this form: 


Ck{t = nT) 


2 _ ^—cnk^ 

Ktk^ 


1 

Kl (fc2 + ’ 


( 20 ) 


where c, Ki and Kt are some constant ^ 0(1). The first and second terms corresponds to transverse and longitudinal 
fluctuations, respectively. For finite time (n), the trasverse term accounts for the linear Floquet divergence of the 
correlation at /c = 0 [Fig. |^b)] in the main text . k^^ = ^ in the longitudinal term is the correlation length which 
diverges near the continuous phase transition [Fig. |^d)]. 

The correlation length near the LC-Ui boundary diverges as ^ ^ This property results in a Ginzburg 

criterion in d-dimension: 


{AS^f W 


«1, 


( 21 ) 


corresponding to a upper critical dimension of 4 in the LC phase. Thus, despite its dynamical properties, the LC 
phase and its phase transition can still be characterized by a divergent correlation length and critical dimensions. 


FROM NORMAL FORM TO THE LOWER CRITICAL DIMENSION 


While the normal-form description (see main text) shows that the limit cycle has a Goldstone mode, and the 
Gaussian analysis above shows that its lower critical dimension is two, the link between these two observations is not 
explicit. Here, we address this gap. This is a standard calculation in the Keldysh formalism; we indicate the main 
steps, and refer the reader to the review article by Kamenev and Levchenko [4] (whose notation we follow) for details. 
We proceed by rewriting the normal-form equations using the Martin-Siggia-Rose approach mu- We then transform 
to the rotating frame, and keep only the phase fluctuations of A (i.e., A(x, t) = Aq exp(i^(x, t))). In terms of these 
phase fluctuations, the quadratic action (written in the classical-quantum basis [4]) is 

Og) = Al [eg{dt - ^ (22) 

where Oc^q is the classical (quantum) phase field and represents the Langevin noise amplitude. From this it is 
straightforward to compute the phase correlation function (^k^-k), as this is the Keldysh component of the 0 Green’s 
function computed from the action above. This gives the result quoted in the main text. 
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PHOTON CORRELATIONS 


To compute the photon power spectrum [Eq. in the main text, we derive the exact equation of motion for the 
two-time correlation function from Eq. 0 : 




= + r)) - r“(cT+(i)cT“(i + r)) - 7<5“^(a+(t)) 

+ 2 E + T)al{t + r)). 

fl 


(23) 


We focus on the phase space where spatial fluctuation is negligible. In this case, the mean-field decoupling of the last 
interaction term (i.e. + t)) ~ + ^))) results in a differential equation 

with time-dependent coefficients. We numerically solve this equation to obtain the power spectrum. The fact that 
the Eq. (23) shares the same structure as Eq. 0 and is usu ally referred as the quantum regression theorem E]. 
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